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Abstract. We investigate the spectrum of a typical non-self- 
adjoint differential operator AD — —dP/dx^ (g) A acting on 
L^(0, 1) (8) C^, where A is a 2 x 2 constant matrix. We impose 
Dirichlet and Neumann boundary conditions in the first and second 
coordinate respectively at both ends of [0, 1] C K. For A e M^^^ 
we explore in detail the connection between the entries of A and 
the spectrum of AD, we find necessary conditions to ensure sim- 
ilarity to a self-adjoint operator and give numerical evidence that 
suggests a non-trivial spectral evolution. 

1. Introduction 

In this paper we investigate spectral properties of the hnear operator 
AD acting on L^(0, 1) ® where A is a 2 x 2 constant matrix and D 
denotes the ordinary differential operator 

\ __f<P"\ 0(0) = 0(1) = 

W) ' 7'(0) = 7'(1) = 0. 

The apparently simple combination of Dirichlet and Neumann 
boundary conditions allows self-adjointness if, an only if, A is real 
and diagonal. If A is non-diagonal and upper-triangular the numer- 
ical range of AD is a large sector of C Otherwise it is the whole of C 
preventing us from applying the theory of sectorial sesquilinear forms 
in a straightforward manner. 

Our main goal is to explore the connection between the entries of 
the matrix A and the location of the spectrum of AD in the complex 
plane. In p| R. F. Streater considers the particular case 



D 



^ = ( 1 T 1 , 7 > 0, 



vl/27 h 

in order to find necessary conditions for the stability of small pertur- 
bations about the stationary solution of certain non-linear system of 
parabolic equations. Streater's system represents a thermodynamical 
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model for hot fluid in one dimension and the localization of the spec- 
trum is achieved by constructing a non-unitary transformation that 
makes AD similar to a non-negative self-adjoint operator, hence the 
spectrum of AD is real and non-negative. This similarity transforma- 
tion does not work for other matrices and a slight modification of the 
entries of A can destroy reality of the spectrum (cf. sections so 
the general case should be attacked by other methods. 

Although this paper mainly concerns A G M^^^, the results of sec- 
tions refer to any complex 2x2 matrix. The core results are to 
be found in section |^ where we present an exhaustive description of 
the spectrum of AD in terms of the entries of A. Among various other 
unexpected conclusions, the following three epitomize the complexity 
of the problem to be considered: 

a) When A is triangular and non-diagonalizable, AD is not similar 
to a self-adjoint operator but the spectrum of AD is real (theorem 

b) The spectrum of AD can be non-real even when both eigenvalues 
of A are positive and equal (theorem pT]) . 

c) There is a continuous family of matrices A whose eigenvalues do 
not intersect the real line but such that the spectrum of AD is 
real (theorem |2^) . 

The last two assertions show that the spectra of A, D and AD are 
typically unrelated. 

The crucial idea in section ^ is to reduce the four-parameter problem 
of localizing the spectrum of AD in terms of the entries of A, to five two- 
parameter cases and describe separately each of these cases. Sections 
are devoted to describing the various properties of AD we will use 
in section ^, whereas section |^ is devoted to numerical computations 
which illustrate some of the results reported. In section Q we find the 
boundary conditions associated to the adjoint of AD and compute the 
numerical range of AD. In section ^ we show that the resolvent of 
AD is compact for all non-singular A. In Section ^ we explore the 
stability of the spectrum of AD in the sense of @] and [H, and provide 
estimates which allow us to enclose the spectrum of AD in angular 
regions when A is subject to various constraints. In section § we use 
standard ODE methods to compute the transcendental function of the 
spectral problem associated to AD. 



2. Definitions and notation 

Let be a linear operator whose domain is denoted by Dom(i^). 
Throughout this paper Spec {K) stands for the spectrum of K and the 
numerical range of K is defined to be 

Num(i^) := {{KfJ) : / G Dom{K), \\f\\ = 1}. 
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We recall that the numerical range of any linear operator is convex and 
that if Spec (X) ^ 0, then 



Spec {K) C Num {!<) . 

If = K* and Spec (A') C (0, cxd), we will say that K is positive 
and write X > 0. \i K ^ K* and Spec(X) C [0, oo), we will say that 
K is non-negative and write > 0. 

Below and elsewhere \v\ denotes the norm of a vector v e C^. The 
norm of any 



/ ^ ^-J G L^(0, 1) ® 
is the standard Hilbert tensor product norm 



(/,/)= f\f{x)\'dx^ [\\<f>{x)\' + \j{x)ndx. 
Jo Jo 



Unless explicitly stated, we denote 

A :-. 



a b 
c d 

The complex numbers a+, a_ denote the eigenvalues of A and the non- 
zero vectors denote the eigenvectors 

Av± = a±v±. 

If a+ and a_ are real and different, we adopt the convention a_ < a+. 
Notice that the v± are not necessarily orthogonal. 

Let W^'^ be the Sobolev space of all / e L'^{0, 1) ® such that 
the generalized derivative /" e L^(0, 1) (8) C^. We define rigorously the 
domain of AD as 

Dom(L>) = {/ e W2'2 : (f)(0) = 0(1) = 0, 7'(0) = = 0}. 

If A is invcrtible, it is standard to show that AD is always a closed 
densely defined hnear operator acting on L^(0, 1) C^. 

Lemma 1. If A is singular, then AD is not closed in the domain 
Dom{D). 

Proof . Let f G be a non vanishing vector such that Av = and 
let f{x) := vx G L2(0, 1) ® C^. Clearly / ^ Dom(D). Let 0„ be a 
sequence of smooth functions whose support is compact in (0, 1) and 
such that (j)n{x) — > x in L^(0, 1). Then 0„v G Dom(D) and (t)nV — > /. 
Also 

AD(t>r,{x)v = -(t>'l,{x)Av = 0, 

so that AD{(j)riv) is a convergent sequence in L^(0, 1)0C^. We complete 
the proof by noticing that if AD was closed, then we would have / G 
Dom(i:>). ■ 



LYONELL S. BOULTON 



For the rest of this section and in sections we will assume without 
further mention that A is non-singular. In section |^ we will consider 
again singular A. 

In order to show that AD is in general non-self-adjoint, let us com- 
pute the adjoint (AD)*. Let 

Then the boundary conditions for D can be rewritten as 

P/(0) = P/(l) = 0, (/ - P)/'(0) = (/ - P)/'(l) = 0. 

Lemma 2. The adjoint of AD is 

{ADYf = -A*f", 
for f G W^'^ subject to the boundary conditions 

P/(0) = P/(l) = 

(/-p)/(o) = (/-p)r(i) = o 

where P = is the rank one projection such that 
Ran (^P) = Ran {A{I - P)) ^ 

Ran [l- P) = Ran(AP) ^. 
Proof . For / G Dom(P)) and g e 1^(0, 1) ® C^, 

{ADf,g) = - [\Af"{x),g{x))dx 



(1) 



(2) 







•1 

|0 



{APf',g)\'i+ / {Anx),g'{x))dx. 







We ought to find a complex 2x2 matrix B and impose boundary 
conditions on g, for 

{f,{ADrg) = - [\fix),Bg"{x))dx 



{B*f,gX+ / {B*nx),g\x))dx 
Jo 

{B*il - P)f,gX + /' {B*f\x),g\x)) dx. 
Jo 



and 

{ADf,g) = {f, [ADYg). 

This must be true in particular for all / and g with compact support 
in (0, 1), so clearly B = A*. 
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Let the boundary conditions for (AD)* be given by (|l|) where P = 
is a non-necessarily orthogonal projection on C^, we show (^. If f,g 
are smooth functions supported in [0, 1), then 

APf'iO), (/ - P)g{0)) = (a{I - P)/(0), Pg'iO) 



where /(O), /'(O), g{0) and g'{0) are arbitrary vectors in C^. If /'(O) 
0, the right hand side should vanish for all f{0),g'{0) G C^, so that 



Ran(^p) = Ran {A{I - P)) 



If /(O) = 0, the left hand side should vanish for all /'(0),^(0) G C^, so 
that 

Ran(^/-P) = Ran (AP) ^. 
Since A is non-singular these two spaces are one dimensional. ■ 

Corollary 3. AD is self-adjoint, if and only if A is real and diagonal. 

Proof . Using the notation of lemma 0, AD is self-adjoint, if and 
only if 

A = A* and P = P. 
The latter occurs, if and only if 

These conditions ensure A real and diagonal. ■ 

We now show that due to the boundary conditions we have chosen, 

Num(AP)) =C 

for a large family of non-diagonal matrices A. This prevents us from 
employing the theory of sectorial sesquilinear forms in order to find the 
spectrum. 

Theorem 4. Let A be a non-singular matrix. 
a) // A is an upper triangular matrix ( that is c = 0), then 



Num (AD) = {rz -.r e [0, oo), z G Num (A) }. 
b) // A is not an upper triangular matrix ( that is c ^ 0), then 

Num (AD) = C. 
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Proof . Since is always an eigenvalue of AD (cf. section then 
e Num [AD) . For / G Dom(L'), 



{ADf, f) 



{Af"{x)J{x))dx 
{Af'J)\U f\Af'ix)J'ix))dx 



i ' V7 



Oy 'V7 



^1 

+ / {Af{x)J\x))dx 

1 JQ 

+ [\Af{x),f'{x))dx. (3) 
1 Jo 



Case a): call 



$ := {rz : r G [0, oo),z E Num (A) }. 
Then $ is a convex set and 



$ = {rz : r G [0, oo), 2; G Num (y4) } 



{{Av,v) : i; G C2}. 



If c = 0, 



so that 



A 



0/'V7, 



j ' U 



{ADfJ)= / (A/'(x),r(x))da;. 
This and the fact that $ is closed and convex, yield 



Num (AD) C $. 

In order to prove the reverse inclusion, let f G be such that |f | = 1 
and let 



z := {Av, v) G Num (A) 



For alH > 5, let 



' ^-7(^*^/2) if 0<x<2/t 
A/4-io/t ~ ~ 

^,(x):=<! ^ if 2/t<x<l-2/t 

l-co.(^t(x-l)/2) l_2/t<X<l. 



V4-10/t 

Then Vt(0) = V-tll) = ^I(O) = ^1(1) = 



[ \^t{x)\'^dx = 1 and /" \'^[{x)\'^dx 
Jo Jo 
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Let ft := vipt £ Dom(Z}). By construction \\ft\\ = 1 and 

{AD ft, ft) = f{Af[{x)J',{x))dx 
Jo 

= {Av,v) [ \^'t{x)\'dx 
Jo 



8t- 20 

Thus by taking t —>■ oo, from the fact that G Num [AD) and since 
Num {AD) is convex, we gather 



Num (AD) D $. 

Case b): now c 7^ 0. Let 2; be a fixed non-zero complex number. Our 
aim is to find functions fs G Dom{D) parameterized by £ > 0, such 
that II /ell = 1 and {Af^, fe) is close to z for small e. 

For < e < 1/2, let 

{^sin(a;7r/e) if 0<x<e/2 

^[l-cos(27r(a;/£-l))] if e/2<x<e 
if e < X < 1. 

Then, straightforward computations show 0e(O) = (pdl) = = 0, 

0^(0) = c-\ 



1 ll^S rl ^ 

|0,(x)|Mx = and \4>'^{x)\Mx = —. 

For all £ > small enough, we define the required test function as 

where 



aie) 



, / llb|2<r3 

:= Vl - I.| W = f - lilir 

is independent of x. By construction G Dom(i5) and 

\\f,r = \\z<P,r + a{ef = l. 

According to @, 

{ADf^J,} = + 



0; ' 10 



za{e) + (Ail],(l]) [' \z\'mx)\'dx 



as 



Since a{e) 1 as e ^ 0, the above shows {ADf^, f^) ^ z as 
£ — i> 0, so that z is an accumulation point of Num {AD) . By moving 
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2 G C, any complex number is accumulation point of Num (AD) . Since 
Num [AD) is convex, the only possibility for Num {AD) is to be the 
whole complex plane. ■ 



3. The resolvent of AD 

In this section we show that the resolvent of AD is compact for all 
non-singular A. In general it is false that the product of a bounded 
operator and an operator whose resolvent is compact has compact resol- 
vent, however if we know in addition that the spectrum of the product 
is not the whole of C, then the assertion is true. 

We first show that the resolvent of D is compact by making use of 
its self-adjointness. Since the constant function 

is in Dom{D) and ADfo vanishes, 

G Spec (AD) . 

Proposition 5. If A is a diagonal matrix, then 

Spec (AD) = {a.Tr^Ai^a+Tr^A;^}^^^. 

The zero eigenvalue is always non- degenerate and all the remaining 
eigenvalues are of multiplicity no greater than 2. 

Proof . Let /o G Dom(D) be as above. For all = 1, 2, . . . , let 

Then fk G Dom(D), 

ADf2n-i = (a+7rV)/2„_i, ADf2n = (a-vr^n^)/2„ 

and {/fcjfcLo is a complete orthonormal set in L^(0, 1) ® C^. ■ 
According to corollary |] and the above proposition, D = D* > 
and 

Spec(D) ={n'e}Zo- 

Since the eigenfunctions {fk}'kLo form a complete orthonormal set, the 
resolvent of D is compact. 

Let us now rule out the possibility Spec (AD) = C. 
Lemma 6. For any non-singular A G C, 

Spec (AD) ^ C. 
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Proof . Fix the matrix A. Since 

AD - X = A{D - XA-^), 
the complex number A G Spec (AD) , if and only if 

G Spec{D - XA-^) . 

Let H{X) := D — XA^^. Then the family of operators H{X) with 
domain Dom(D) independent of A is a holomorphic family of type 
(A) for all A G C. Since is a non-degenerate isolated eigenvalue of 
H{0) = D and A~^ is bounded, there exist an open neighbourhood 
G t/ C C such that H{X) has a non-degenerate isolated eigenvalue, 
(denoted by /i(A)) close to for all A G t/ and /i(A) is a complex valued 
holomorphic function in U (cf. []^, th.XII.8]). 

If there exists some Aq G f/ satisfying ji^Xo) ^ 0, then ^ Spec (if (Aq)) 
so that Aq ^ Spec (AD) . Hence, in order to show that Spec (AD) ^ C, 
it is enough to show that /i ^ 0. For this we find the first coefficients 
in the Rayleigh-Schrodinger series expansion of /x about 0. Let 

/i(A) = /io + yUiA + /i2A^ + . . . X eU. 

Since /i(0) = 0, = 0. Since ||/o|| = 1 and H{0)fQ = D/q = 0, we 
compute directly /xi (cf. P, remark 2.2, p. 80]) by 

Ail = (A-^foJo) 

fo^ 



A- 



1 ' \1 



If A is such that a ^ 0, 



so that fii does not vanish and hence yU ^ 0. 

Let A be such that a = 0. Then /xi = so we compute fi2. Let fk be 
the eigenfunctions of D as in (|) so that = 1 for all A; = 1, 2, ... . 
Let A2n-i = A2n := TT^n^ for all n = 1, 2, . . . so that 

H{0)h = Dfk = Xkfk. 
Then (cf [|, remark 2.2, p.80]) 

{A'^foJk) {A-'hJo) 



-/i2 



fc=i 



We compute each term in the series. Since a = and A is invertible, 
then b and c do not vanish and 

,_(-d/ibc) l/c\ 

^ -[ 1/b o)- 

Hence 

{A-'fo, fk) = (a-' (J) , fk{x)^ = ( Co"") ' ^^(^)) 
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SO that 

• 1 



and 



(A-V../.„)^V2/(f/').(,„,(L)))d. = 



2/c / sm{'n-nx)dx 
Jo 

if n = 2m 

2V2/{cTTn) if n = 2m-l 

for m integer and n — 1,2, ... . On the other hand 



SO that 



and 



\/2/{b) i sin(7rna;)da; 



if n = 2m 

2%/2/{lmn) if n = 2m - 1 

for m integer and n — 1,2, This yields 

{A-'fo, fk) {A-'fk, fo) = I 8/(^c7rV) if kttZ-l 
for m^ 1,2,.... Thus 

"^^ = ^ 6c7r4(2m - 1)^ ^ ° 

m=l ^ ' 

SO that // ^ as we required. 



Theorem 7. For all z ^ Spec (>1-D) , the resolvent {AD — z) ^ is com- 
pact. 
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Proof . Since D is non-negative and it has compact resolvent, 

AD + A = A{D + 1) 

has a compact inverse. Let z Spec (AD) , then 

AD-z = AD + A-A-z 

= {l-{A + z){AD + A)-^){AD + A). 

Hence 

{AD + A)-^ = {AD - zy\l -{A + z){AD + A)-^), 

so that 

{AD-zy^ = {AD-zy\A + z){AD + Ay^ + {AD + A)~^ 
= {{AD - z)-\A + z) + 1){AD + A)-\ 

Thus {AD — z)~^ is compact as needed. ■ 
Theorem |^ shows that the spectrum of AD consists entirely of iso- 
lated eigenvalues of finite multiplicity. Since the eigenvalue problem 
ADf = A/ is a constant coefficient system of second order ordinary dif- 
ferential equations, due to the fact that we have a combination Dirich- 
let and Neumann boundary condition at both ends of the interval, the 
multiplicity of each eigenvalue is never greater than 2. 



4. ASYMPTOTICS OF THE RESOLVENT 

We now investigate the asymptotic behaviour of the resolvent norm 
of AD. The results we discuss in this section are connected with the 
stability of the heat semigroup e~^^*. They are also relevant from the 
computational point of view and they are closely related to both local 
and global stability of the spectrum (cf. [^, [^] and the reference 
therein). The present approach is motivated by analogous reports on 
non- self- adjoint Schrodinger operators in 0, and [^]. 

Let 



J :-- 



1 
-1 



Below and elsewhere we will denote hj D := JD. According to lemma 
D = D*. According to lemma ^, 



Spec(D) ={±7rV}- 







each eigenvalue being of multiplicity 1. We will employ part b) of the 
following theorem in the proof of theorem [21|-b). 

Theorem 8. Assume that there exists a non-singular diagonal matrix 
B such that R-^AB = {B'^AB)* > 0. Then 

a) AD is similar to a non-negative self-adjoint operator. 
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b) AD is similar to a self-adjoint operator whose numerical range is 
the whole real line. 

Proof . Let C := B'^AB so that C = C* > 0. Since diagonal 
matrices commute with the boundary conditions, AD is similar to CD. 
For the same reason and since diagonal matrices also commute with J, 
AD is similar to CD. 

By hypothesis, the square root C^^^ = (C^^^)* > 0. Then 

CD = C^I\C^I''DC^I'')C-^I'' = C^I'^KC-^I'', 

where 

K = C'/^DC'/' 

Bom{K) = {f e L\0, 1) ® : C^^^f G Dom{D)}, 

so that CD is similar to K. Since D = D* > 0, then K = K* > 0. 
Analogously CD is similar to 

k := C'^'^DC'^' 

where Dom(i^) = Dom(fC). Since D = D*, then K = K* . Further- 



Num [D 

and 

'kfj) = (DC'/'f,C'/'f 



the numerical range of K is the whole real line. ■ 
Let A be as in the hypothesis. The similarity to a self-adjoint oper- 
ator ensures the existence of a constant fc^ > 1 such that 

II (^^ - < z Spec (AD) 

dist [z, [U, oo)) 

and 

\\(AD-z)-^\\<—^\— z^Spec(AD 
^ II - dist(z,M) ^ ^ y 

These identities show that although the numerical range of AD and 
AD are in general the whole complex plane, the eigenvalues of these 
operators are stable in the sense of 0. 

If we assume the weaker condition C + C* > 0, we show how to 
recover part of the above estimate. We start with a preliminary lemma. 

Lemma 9. Let A be such that Num (A) C {Re(2) > 0}. Then 
Spec (AD) C {Re{z) > 0} and there exists k > independent of z, 
such that 

\\{AD - z)-'\\ < -1^ Re(^)<0. (5) 



' rz). 
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Proof . Let r > and let z ^ [0, oo). Then 

AD - z =A{D - zA-^) 

=A [{D - rz) + z{r - A' 

=A [l + (r - A-'^)z{D - rzy^] {D 

Therefore z ^ Spec (AD) , whenever 

\\{r-A~^)z{D-rz)-^\\<l. (6) 

We show that there is always r > independent of z, such that this 
holds for all Re{z) < 0. 

Since D > and e Spec {AD) , 

IIP--rNl = ;:^- 

Thus 

||(r - A-^)z{D - rz)-^\\ < ||1 - r-^A-^\\. 
The hypothesis we imposed on A is equivalent to saying 

A + A*>0, 

then 

A-^ + {A'^)* = A~\A* + A){A'^)* > 0. 

For all V eC^, 

-r'^A-')vf = {{I -r-\A-' + {A-^)*)+r^\A-^)*A-')v,v) 

= \v\^- r^^ { {A-^ + (A-^)* + r'\A-^)*A'^) v, v) 

Hence there exists a constant fco > independent of r (and z), such 
that 

\\I -r-^A-^\\ < l~r'^ko 

when r is large enough. For such an r, identity (|^) holds for any 
Iie{z) < 0. This shows that Spec (AD) must be enclosed in the right 
hand plane. Furthermore 

\\iAD-z)~'\\ < \\A-'\\\\{D - zA-Y'\\ 

< \\A-^\\\\{D-rz)-'\\ {l + z{D-rz)-\r-A-')y' 
II 4-111 

< }L^J2\\z{D -rz)-\r - A-^)\\^ 

^^^^ 1=0 



< MZ^y 111 _ ^-1^-1 



1=0 



k 

< — 

\z\ 



so d^) is also proven. 
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Below and elsewhere we denote by fl the set of non-singular diagonal 
matrices and 

S{a, (3) := {z eC : a < arg(2) < (3} a < (3. 

Theorem 10. // there exists B eVL such that 

^Mm [B-^ AB) d S {a, 13) (3-a<7T, 

then Spec (AD) C S{a,l3) and for any small enough e > there exists 
fcg > independent of z, such that 

\\{AD - z)-^\\ < z^S{a-e,p + e). 

Proof . Let C := B'^AB, so that 

Num(C) c S{a,(3). 

Since B commutes with the boundary conditions, AD is similar to 
CD and so it is enough to show the theorem for CD. Now, for all 
-(a + 7r/2) <^<7i/2-P 

Num (e^'^C) C {Re{z) > 0}, 

so we just have to apply lemma |^ to e^'^C. ■ 
The constant fee of this theorem is in general strictly greater than 
1, therefore this is weaker than the similar condition for m-sectorial 
operators in p. 279]. 

If A is triangular, the hypothesis of the above theorem does not 
necessarily hold. For instance if 

then 

Num (A) = {a + z:\z\< 1/2} 

and so for small a the numerical range contains the origin. Nonetheless 
by using a similarity transformation and an approximation argument, 
we can show positivity of the spectrum whenever both of the eigenval- 
ues of A are positive (a > in our example). The conclusion about the 
spectrum of the following result will be improved in theorem |20|. 

Corollary 11. Let A be either upper or lower triangular. Ifa>d>0, 
then 

Spec (AD) C [0,oo) 
and for all e > there exists fc^ > independent of z, such that 

\\{AD-zn\<^^ 

for all z ^ S{—e,e). 
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Proof . If A is upper triangular the proof is similar so let us assume 
that 



A 



a 

c d 



Let 

rj^\0 r^^J \rc d 
Then AD is similar to A{r)D for all r 7^ 0. Put 

C(r) := A{r)+A{r)* 

Then C(r) = C(r)*. The eigenvalues of C{r) are 



2a rc 
rc 2d 



a + d± v/(a-c/)2 + r2|c|2, 

thus for small r > 0, C(r) > 0. The numerical range of A{r) is an 
ellipse with focus at a, d and principal axis in the vertical direction of 
the order of r. By taking r —>■ 0, theorem ^ completes the proof. ■ 

If A is as in the hypothesis of corollary O, there does not exist 
Ben such that B-^AB = {B'^AB)* or B^AJ)B = {B-\AJ)B)* 
so theorem |^ is not applicable. We show that at least in one case AD 
fails to be similar to self-adjoint. 



Theorem 12. Let 



A = { , 1 a > 0. 

1 a 



Let e > and z{r) := Aan'^r^ ± ie. Then there exists a constant > 
independent of r, such that 

\\{AD - z{r))-^\\ > ky/^ r = 1,2,... . 

Proof . Fix e > and let z(r) := Aan'^r'^ — ie. Without loss of 
generality we can assume r = 3,4,.... Throughout the proof the 
constants Ij are assumed to be positive, possibly depending upon e but 
independent of r. In order to show the desired conclusion, it is enough 
to find G Dom(D) and /q, such that 



\\ADfr- Z{r)fr\\ ^ , _i/2 



II/. 



< kr-'" (7) 



for all large enough r. 
Let 



f=\Z^] GDom(Z}). 
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Then 



- ^ - (t (^::) - .(D 



7 



-ad)" — zir) 



III 



07" — z{r)'-) j 



Hence 



2 _|_ ll,,l|2 



and 

\\ADf - z{r)ff = ||a0" + ^(r)0f + \\a^" + ^(r)7 + ^''f- 
We now define the appropriate fr G Dom(_D) satisfying (|^). Let 

7^(0;) := cos(27rrx) 

Then ||7^f = 1/2. Let 

—ie cos{2Txrx) / {Aix'^r'^) if a; G (1/r, 1 — 1/r) 
\ if x ^ (l/(2r), 1 - l/(2r)) 

be such that 0^ is smooth and 

a) \(t)r{x)\ < £/(47rV2) for all a; G [0, 1], 

b) WXx)] < h/r for all x ^ (1/r, 1 - 1/r), 

c) < /a for all x ^ (1/r, 1 - 1/r). 

Then 

||0^f = f |0,(x)|2dx < £V(167rV) < hr-" 
Jo 

and 
Hence 



/. = C^^^ ) G Dom(D) 



and 



l/2<||Mr=||0.f +||7r.f <1 

for all large enough r. Ifl/r<x<l — 1/r, 

a(j)"{x) + z{r)(j){x) = a(f){x)" + 4a7r^r^0(a;) — iecpi^x) 

= [aie cos(27rrx)" + Aan'^r'^ie cos(27rrx) + 

— e'^ cos(27rrx)]/(47r^r^ 
= — cos(27rrx)/(47r^r^). 
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Then, a) and c) yield 

\\af + z{r)(f)f = [ \a^"{x) + z{r)(f){x)\^dx 
Jo 

< / k/r^dx+ / k + lr/r^dx 

< kr-^ + hr-^ + hr-\ (9) 

Also, 

07" (x) + z(r)7(x) + = a'j"{x) + 4a7r^r^7(x) — ie'y{x) + 

= a cos{2TTrx)" + 4a7r^r^ cos(27rra;) + 

— ie cos(27rrx) + (f>"{x) 
= (p" (x) — ie cos{2Tcrx). 

Then for 1/r < x < 1 — 1/r, 

a7"(x) + z{r)'y{x) + (f)"{x) = ie cos{2nrx) — ie cos{27Trx) = 
and thus c) yields 

fi 







|a7" + ;z(r)7 + 0"||^ = / \ai' (x) + z{r)-f{x) + (f)" {x)\Mx 

\ie cos{2Trrx) — (f)"{x)\'^dx 



x^[l/r,l-l/r] 



Jx^[l/r,l-l/r] 

< kr-\ (10) 

In order to complete the proof for z{r) := 4a7rr^ — ie, notice that 
(D, (D and (|10|), show (0). On the other hand, if z{r) := Aairr'^ + ie it 
is enough to substitute (j)r by — 0^ and repeat the above computations. 

■ 

This result is still valid for 

'a 1 



,0 a I 

Indeed, it is enough to put := sin(27rrx), 

±iesin(27rrx)/(47r^r^) if x G (1/r, 1 — 1/r) 



^^^^y-'- ^ if (l/(2r),l-l/(2r)) 

and carry out similar calculations. Since the resolvent norm of self- 
adjoint operators remains bounded in horizontal lines, the above AD 
can not be similar to any self-adjoint operator. 

Let Qr be the set of all non-degenerate real diagonal matrices. If A 
does not satisfy the hypothesis of theorem (for instance the numer- 
ical range of A is an ellipse centered at the origin), but A is "close" in 
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some sense to an alternative to theorem can be established. We 



will employ this result in the proof of theorem |20 . 

Theorem 13. Let there exist B E Qr such that 

\\AB-I\\ < 1. 

Let uj := arcsin(||y4i? — /||) with < uj < 7r/2. Then 

Spec {AD) C S{—uj, uj) U 5'(— 7r — cu, — tt) 

and for any small enough £ > there exist fc^ > independent of z, 
such that 

\\iAD-z)-'\\<^^ 

for all z ^ S{—uj ~ e,uj + e) U S'(— tt — uj — e,uj — n + e). 

Proof . If = 0, y4 G fir so the conclusion is consequence of corollary 
|. Let UJ > 0, let / := \\AB - I\\ and put C := B'^ G Q,-- Then 
CD = (CD)* and according to the hypothesis < / < 1. 

Let 2; G C be such that z ^ S{—uj, uj) U S{—7[ — uj,uj — tt). Then 

{AD-z) = CD + {A-C)D-z 

= [I+{AB-I)CD{CD-z)-\CD-z). 
Since CD is self-adjoint and by definition w = arcsin(/), 

\\{AB - I)CD{CD - z)-^\\ < l\\CD{CD-z)'^\\ 



< I sup 



X 



X — z 



< sup^-^<l, (11) 
xm 1 1 — I 

so that 

[/ + [AB - I)CD{CD - z)-^] 

is invertible. Hence 

z ^ Spec (AD) 

and 

{AD - z)-^ = {CD - z)-^[I + {AB - I)CD{CD - z)-^]-^ (12) 

for all < I arg(2;)| < tt. This encloses Spec (AD) . 
In order to show the second part, let 

z ^ S{—uj — e,uj + e) U S'(— tt — uj — e,uj — n + e), 

for small e > 0. Then there exist a constant li{e) > independent of 
z, such that 

\\{CD-z)-'\\<^$. 
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Also, there exist a constant < 12(e) < 1 independent of z, such that 

^ < hie). 



sup^ T, 

xm J- — - 



These two estimates, (|ri|) and (|I2D yield 



WiAD - z)-'\\ < 



This shows that if A„ G C^^^ is a sequence of non-singular matrices 
and there exists B &VLr such that 



as n — >■ 00, then 



\An-B\ 



Spec(y4„i:') 







Corollary 14. Let A he either upper or lower triangular. If a,d & M. 
and ad < 0, then 

Spec (AD) C M 
and for all e > there exists > 0, such that 

\\{AD - z)-^ < 

for all z ^ S{—e, e) U S'(— vr — e,e — n). 

Proof . It is similar to the proof of corollary |11[ Assume without 
loss of generality that 6 = and let 



A{r) 



1 OA A 

r \0 



a (T 
rc dj ' 



Then AD is similar to A{r)D for all r > 0. Put 





then 



C 



P(r)C-/| 



d- 



( ° 




\rc/a 





r\c/a\. 



Let uJr '■= arcsin(r|c/a|). According to theorem [T3|, for all < r < \a/c\ 
Spec (AD) = Spec{A{r)D) 

C S{—UJr, UJr) U S{ — Tr — UJr, UJr — n) . 

By taking r small enough, theorem |13| yields the desired estimate for 
the resolvent norm. By taking r —>■ 0, a fortiori Spec (AD) C M. ■ 
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5. The Hamiltonian ODE system 

In this section we find an entire function wliose zeros coincide witli 
Spec (AD) . Tfiis is made by computing tlie transcendental function 
associated witfi tlie 2x2 system of ordinary differential equations as- 
sociated to AD via standard ODE arguments. 

Let the 2x2 constant coefficients second order eigenvalue problem 

-Af" = X'f (13) 
P/(0) + (/ - P)/'(0) = 

P/(l) + (/-P)/(l) = 0. ^ ^ 

We will say that the complex number A is an eigenvalue of the system 
(|T3[)-(0), if there exist a non- vanishing / G C°°(0, 1) ® satisfying 
(|T3p and the boundary conditions (|T^). By regularity, is an eigen- 
value of AD, if and only if A is an eigenvalue of (pTsD -fp^. Our aim is 
to find a holomorphic function, denoted by EV{x) below, whose zeros 
coincide with the eigenvalues of ([I3|)- ([I^) . 

We proceed in the classical manner. Let the decomposition in Jordan 
canonical form of A be 

A =: VCV-^ 
where the Jordan matrix C is either 

C=h M or C=h ° 
U a- J VI a+ 



and 

V :-- 



Vi V2 



Then (|T^) - ([l^) is equivalent to the 2x2 system 

-Cg" = X'g (15) 
PVg{0) + (/ - P)Vg'iO) = 
PVg{l) + {I-P)Vg'{l) = 0. 



(16) 



In order to solve (p!5|)-([T6|), we reduce it to a first order 4x4 system as 
follows. For all A G C, let 



and let 



By regarding 



Pa = ( Jc-i o) e C'^' 



V3 V^j 



$= (|,) GC^ 
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one sees that ([l5|)-(p!6D is equivalent to 

$' = Bx<^ (17) 
^$(0) = ^$(1) = 0. (18) 

In order to solve (|l^-(|18D we must find a fundamental system of 
solutions. Let 61,62,63,64 be the standard orthonormal basis of the 
Euclidean space C^. A straightforward computation show that 



exp{Bxx)ej 



xe[o,i], J = 1,2,3,4 



is indeed a linearly independent fundamental system for (|T^)-(|1^). 
Hence, A is an eigenvalue of this system, if and only if there exist 
fci, ^2, ^3, ^4, such that 



$(x) = kj exi)(Bxx)ej 



(19) 



is non-vanishing and satisfies the boundary conditions. 

We now proceed to compute EV{x). The exponential of Bxx is given 

by 



exp{Bxx) 



-XC-'/^ sin(AC-i/2x) cos(AC-i/2a;) 



for X G [0, 1]. In theorems |T5| and |T6| below, we split our computation 
into two cases depending upon the Jordan matrix C. 



Theorem 15. When 



C 



a_|_ 
a_ 



A is an eigenvalue of the system (fT^j-dl^, if and only if EV{X) = 
for 



EV{x) := l^Yl 



1 — cos 



X 



COS 



.2„,2 



2„,2 " 



sm 



+ 



sm 



X 



Proof . Notice that EV{0) = 0. Assume A 7^ 0. According to the 
hypothesis. 



^1/2 



1/2 n 

a / 



a 



1/2 



and 



^-1/2 



-1/2 
+ 

a] 





-1/2 
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COS 



Xx 



/a+ . Xt 

sm 



COS 



Xx 



\ 

Xx 



sm 



A ^ , Ax 



COS 



Ax 







V 











COS 



Xx 



J 



Let be a particular solution given as in (|19]), where the complex 
parameters kj are to be determined. Then 



^$(0) 



kiVi + k2V2 
k3V3 + k4V4 



and 



Ai 



^$(1) 



^ ' k2V2 COS ^ 



t;i cos ^= + ft2V2 COS + 



sm 



A 

A „, A 



—kiv^-^ sin ^= — ^2^4 



sm 



+ k4V2 
X I 



A 



sm ■ 



v 



+ /i;3t>3 COS + k^Vi COS y 



The solution $ satisfies the boundary conditions (|T^), if and only if 

fciVi + fc2f 2 = 

ksVs + kiVi = 
kiVi cos + ^2^ 2 cos ■ 



sm 



A:, 



1)2 wa- 



sm 



-/ti^sm^ 



sm 



+ ^3^3 COS ^= + A;4t>4 COS 



0. 



The determinant of this 4x4 system of linear equations in kj is precisely 
EV{\). ■ 



a+ 
1 



Theorem 16. When 

C 



A is an eigenvalue of the system (fnj-diqj, z/ anc? on/?/ if EV{X) = 
/or 

:=('-l4-)x^-( det V + 



4ai /V 2a 



Proof . Notice that EV{0) = 0. Assume A 7^ 0. One can verify 
directly that 



2^/0:;: + / \ 1 2a 
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Then the four 2x2 blocks of the matrix exp{Bxx) are 

cos^ 



cos(AC-^/2x) = xx 



2a 



3/2' 



/a+ 

sin 



Xx 



COS 



Xx 



sm{XC-^/^x) equal to 



v 



1/2 



sm 



Ax 



: sm 



Xx 



and sm{XC-'^^'^x) equal to 



cos 



1/2 




■ sm 



-m-sin — 



A 



^^TTjv sin + ^ cos 



Let be a particular solution given as in (p!9D, where the complex 
parameters kj are to be determined. Then 

^kiVi + k2V2 



^$(0) 



and 



where 



^^^^^ U$(l)2 



*$(l)i = /ci (t;icos^ + -^sin ^ 



+A;3(^sin^ + 



i'2 



A 

and 

^$(1)2 = fci 
-fcs^sm 



'a+ 2A^/a+ 



sm 



/oif: 

A 



+ k2V2 COS 



/a+ 2a+ 



^COS^ 



i4 — 



sm 



^ sin ^ + ^ sin ^ + ^ cos ^ 



2a 



2a4 



+ 



^ : + A;3 f ^3 cos + sm ^ 



2a , 



/i:4f 4 COS ■ 



The solution $ satisfies the boundary conditions (|T^) , if and only if 

/Cif 1 + k2V2 = 

/jest's + k^Vi = 
^$(1)1 = 
^$(1)2 = 

A rather long but straightforward computation shows that the deter- 
minant of this 4x4 system of linear equations in kj is EV{X). ■ 

We show that AD can have non-real eigenvalues even when the spec- 
trum of A is positive. 
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Example 17. Put 

A : 



( 2/5 + 3i/10 3/5-3i/10 \ 
\^3/20 + 3i/10 17/20 - 3^/10 J " 



Then the eigenvalues of A are a+ = 1, a_ = 1/4, and the eigenvectors 



1). «--a% 



Thus 

= 4i(l - cos(a;) cos(2a;)) = 4i(l - 2 cos^(a;) + cos(a;)) 
so that EV{X)=G, if and only if 

cos(A) = 1 or cos(A) = -1/2 ± i/2. 

Hence 

Sy>cc{AD) = {4fc27r2, (A± + 2kT^f}k& 
where A± = arccos(-l/2 ± z/2) « 2.02 ± 0.53i. 



6. Real matrices 

In this section we explore some connections between the entries of 
the matrix A and the global behaviour of Spec (AD) when A e M^^^. 
Alongside we discuss conditions to ensure similarity to a self-adjoint 
operator. For completeness of the picture, below and elsewhere we 
allow det(A) = 0. 

Our first task is to reduce to two parameters the four that are initially 
given as entries of A. This leads us to five different types of matrices 
to deal with. For a, d e R, let 



a 0\ ■ (a V\ . fa 

' ■= [i d) ' ■= [l d 

a 1\ , . fa —1 

d ^^^=1 d 



We show that the AjD generate any AD, A e M^^^ via similarity trans- 
formations. 



Lemma 18. If A G M^^^, then AD is similar to a AjD for some 



px2 

a,a,d eM. and j = 0, . . . ,4. 
Proof. Let 



A 



If be — 0, the proof is trivial. Let 



a b 
c d 
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aAi 



'he d 

for a = Vbc, a = d/Vbc and d = d/^/bc. If b/c < 0, 

for a = V—bc, a = ±d/^/—bc and d = ±d/ ^ —be. ■ 
The case j = was already described in corollary ^ Indeed if arf 7^ 
then AqD is similar to a self-adjoint operator and 

Spec(Ao£') = {an^-n^.dn^Ti'^}'^^^ C M. 



6.1. Matrix A^. 

Since a and d are real, Ai = A*. Let b± be the eigenvalues of A\. Then 

a + (i± A/(a - + 4 
^ 2 

so that 

i) b^ > b_ > 0, if and only ii ad > 1 and a, > 0, 

ii) b^ <bj^ < 0, if and only if arf > 1 and a, (i < 0, 

iii) 6+ and 6_ have opposite signs, if and only if ad < 1. 

Theorem 19. The following statements are true. 

a) If ad = 1 then Spec(y4iD) = C. 

b) If ad > 1 and a,d > then AiD is similar to a non-negative 
operator so that Spec{AiD) C [0, 00). 

c) If ad > 1 and a,d < then —AiD is similar to a non-negative 
self-adjoint operator so that Spec{AiD) C (— cxd,0]. 

d) If ad < I then Spec(v4iD) C M. 

Proof . If arf = 1, the matrix is singular so according to lemma 
m, AiD is not a closed operator. This shows a). Statement b) is 
consequence of i) and theorem p, and statement c) is consequence of 
ii) and theorem 

Let us show d). For e G M, let 

B{e) := Ai +ie. 

Then 

Num {B{e)) C {\m{z) > 0} e > 

and 



Num(B(£)) C {Im(^) < 0} 



e < 0. 
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According to theorem p!0| , 

Spec{B{e)D) C {Im(^) > 0} e>0 

and 

Spec{B{e)D) C {lm{z) < 0} e < 0. 

Since B{e)D is a holomorphic family of type (A) in a neighbourhood 
of e = and 5(0) = Ai, a fortiori 

Spec(AiD) C M. ■ 

Although Ai = A^, it is unclear to us whether AiD is similar to 
self-adjoint in the latter case. 



6.2. Matrices A2 and ^3. 

Since the results for the matrix ^3 are analogous and shown in a similar 
manner as for A2, we will only consider the latter. 

Theorem 20. The following statements are true. 

a) If ad = then Spec(A2-D) = C. 

b) lfad=^0 then Spec{A2D) = {an'^n'^,dn^n^}^^Q. 

c) If ad > 0, for all e > there exists > independent of z, such 
that 

\\(A2D - z)-'\\ < z^±S{-e,e), 

where the symbol ± is chosen according to the symbol of a. 

d) // ad < 0, then for all e > there exists k^ > independent of z, 
such that 

\\iA2D - zr^ < 



\z\ 



for all z ^ S{—e, e) U S{—7t — e, — tt + e). 
e) If a = d ^ 0, let e > and z^ = Aan^r^ ± ie. Then there exists a 
constant k^ > independent of r, such that 

\\iA2D-Zr)-'\\>K\Zr\'^^ 

for all r = 1,2, ... . 

Proof . If ad = 0, the matrix A2 is singular so according to lemma 
|l|, A2D is not a closed operator. This shows a). 

Let us show b). If a ^ d, the matrix A2 is diagonalizable and 

_ fa-d 0\ fa 0\ f {a- d^^ 0^ 
- ' 1 1 [0 d i-(a-d)-i 1, 



Then, according to theorem 15, 



EV(x) = kn sin ( | sin ( ) 
Vv^/ \VdJ 
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where ko is constant in x. If a = d, A2 is already in Jordan form and 



according to theorem 16 



EV{x) = - 



sm 



X 



Hence in both cases 

Spec{A2D) = {a7rV,d7rV}^^o. 

Statements c) is consequence of corollary ^ and statement d) is 
consequence of corollary For statement e) use theorem 12 and the 
fact that \zr\ is of order r^. ■ 



6.3. Matrix A4. 

Formally speaking, so far the spectrum of AjD for j = 0, . . . , 3 re- 
produces the spectrum of Aj in the following sense: if Aj is non- 
degenerated and both eigenvalues of Aj are positive (negative) then 
Spec (AjD) is non-negative (non-positive), and if the eigenvalues are 
of opposite sign then AjD possess both positive and negative spectrum. 
There is no reason to expect the same for j = 4, in fact this case is 
less simple due to the way the entries of A4 interact with the boundary 
conditions. 

The eigenvalues of A4 are given by 



J± 



(20) 



2. In this case Aa is not a 



iii) 
iv) 



Then 

i) 6+ = b_, if and only if \a — d\ 
diagonalizable matrix. 

b± are real and have opposite signs, if and only if arf < — 1. 

> b_ > 0, if and only if ad > —1, \a — d\ > 2 and a + d > 0. 
6_ < 6+ < 0, if and only if ad > —1, \a — d\ > 2 and a + d < 0. 

v) b± are non-real with b+ = b-, if and only if |a — (i| < 2. 

vi) is singular, if and only if ad = —1. 

Motivated by this and for simplicity, we can divide the plane into 6 
disjoint regions Rk, 

\a — d\ = 2,a ±1}, 
ad < —1}, 



Ri 
R2 
R3 

i?4 

R5 
Rr 



{{a, d) G 
{(a, d) e 
{{a, d) e 
{(a, d) E 
{{a, d) e 
{{a, d) E 



ad > -l,\a - d\ > 2,a + d > 0}, 
ad > —1, |a — (i| > 2,a + d < 0}, 
\a- d\ < 2}, 
ad = —1}. 

Clearly = [jRk- Below we establish the spectral results for A4D 
separately in each region Rk- 
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Two cases are similar to what we have found so far. 

Theorem 21. The following statements are true. 

a) If {a,d) G Rq, then Spec^A^D) = C. 

b) // {a,d) G R2, then Spec{A4^D) C M and A4D is similar to a 
self-adjoint operator whose numerical range is the whole real line. 

Proof . If ad = —1, the matrix A4 is singular so according to lemma 
|I|, A4^D is not a closed operator. This shows a). 
Let us show b). Let J be as in section ^. Then 

A^D = {AiJ){JD) = }^ D = Ad. 

Here A = A* and the eigenvalues of A are 

~ a-d±J(a + dy + 4: 
= 2 • 

Since ad < —1, 6-1- are either both positive or both negative. If they 
are both positive, A > so that theorem ||-b) provides the desired 



SPECTRAL BEHAVIOUR OF A SIMPLE N-S-A OPERATOR 



29 



conclusion. If they are both negative apply the above argument to 

-A^D. m 

In order to find Spec{A4^D) in R/. for k = 1,3,4,5, we ought to 
rely on properties of the transcendental function EV{x). Nonetheless, 
theorem ^ provides some indication of what we should expect, it bases 
on the observation that if both a and d are positive, 

A, + Al= ^ > 0, 



2dj 

so by virtue of lemma H, Spec^A^D) C {Re{z) > 0}. 
Theorem 22. // both a and d are positive, then 

Spec{A4^D) C S{—uj,uj) 
where sinu = 1/^/ad + 1 for < u < n/2. 

Proof . The numerical range of A4 is an ellipse whose foci are 6-1- and 
largest diameter is of length |a — It is easy to see that S{—uj,u!) is 
the minimal sector that contains such an ellipse. Use theorem |10| to 
complete the proof. ■ 

Since 

-1 0\ /-a 0\__/a -1 

y \^ 1 -d) \ o 1) ~ [1 d 

and because of diagonal matrices commute with the boundary condi- 
tions. Spec {A4D) C —S{—LJ, u) where both a and d are negative. This 
also shows that the spectral results for A^D are symmetric with respect 
to the transformation {a,d) i-^ {—a,—d). Below we will employ this 
symmetry often without mention. 

In order to describe Spec{A4^D) in R^, we will make use of the fol- 
lowing technical result. 

Lemma 23. Let a E C be such that Re(a^) > 0, let —1 < c < 1 and 
let 

F{x) := 1 — cos(aa;) cos(aa;) — csin(ax) sin(ax) x G C. 

Then F{x) has an infinite number of zeros in the complex plane and 

a) if c = —1, then F{x) = 0, if and only if sin (Re(a) x) = 0, 

b) if c = 1, then F{x) = 0, if and only if sinh (Im(tt)x) = 0, 

c) if —1 < c < 1, then F{x) only has a finite number of zeros lying 
on the real and imaginary axis. 

Proof . Let a =: p + ifi so that p > p > and let x =: Xi + 1X2 for 
Xi, X2 E M. 
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In order to show a), assume c = — 1. Then 

|F(a;)|^ = \1 — cos[a{xi + ix2)]cos[a{xi + ix2)] + 

+ sm[Q;(a;i + ix2)] sin[a(a;i + ix2) 

— 4[cos^(pa;i) — cosh^(pa;2)]^. 

Hence 

F{x) = 0, 

if and only if cosh(pa;2) = 1 and cos(pa;i) = ±1. This gives a). 
Similarly for b), assume c — 1. Then 

= \l — cos[a{xi + ix2)]cos[a{xi + ix2)] + 

— sin[Q;(a;i + ^0:2)] sin[a(a;i + ix2) 

— 4[cosh^(/xxi) — cos^(/ia;2)]^. 

Hence 

F{x) = 0, 

if and only if cosh(/xa;i) = 1 and cos{iiX2) = ±1. 
Let us show assertion c). If x e R, then 



F{x) — 1 — cos(aa;)cos(aa;) — csin(Q;a;)sin(Q;x) 
= 1 — I cos(«a;)p — c| sin(Q;a;)p 
= 1 — cos^(pa;) — csin^(pa;) — (1 + c) sinh^ [iix] 

and 

F{ix) = 1 — cos{—iax) cos{iax) — csin{—iax) sm[iax) 
= 1 — cos(mx) cos(ma;) + csm{iax) sin(iaa;) 
= 1 — I cos(ma;)|^ + c| sin(i«a;)|^ 
= 1 — cos^dix) + csin^(/xa;) — (1 — c) sinh^(pa;). 
Hence, if — 1 < c < 1, 

lim F[x) = —00 and lim F{ix) = —00. 

»±oo x—*±oo 

Since F{x) is a smooth function, c) follows. 

Finally let us show that F{x) has a infinite number of zeros. Suppose 
that F only has a finite number of zeros 0,Zi,... ,Zn where the Zj 
repeats as many times as its order. Then 

F(x) 



G{x) 



is an entire function with no zeros. By virtue of the Weierstrass factor- 
ization theorem, there is an entire function g{x) such that G{x) — e^^^\ 
Then 



F{x) 



x^ 
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Since it is a combination of sines and cosines, the order (cf. 0, p. 285]) 
in the sense of entire functions of F{x) is A = 1. Thus by virtue of 
Hadamard's factorization theorem, g{x) is a polynomial of degree 1 in 
X and so 

F{x) = p(x)e'=^+' 

for suitable k,l G C. Since p{x) is a polynomial, this is clearly a 
contradiction, so F{x) must have an infinite number of zeros. ■ 

Theorem 24. Let {a,d) G -R5. 

a) // (a, d) e {a^ ~ ad - 1 = 0} n {-2 < a - d <0}, then 

SpeciA^D) = \-k\y[lm{bl^^^)]^} C (-00, 0]. 

b) If {a,d) e {a^ - ad - 1 = 0} n {0 < a - d < 2}, then 

Spec{A^D) = |A;V/[Re(6;'/')]2| C [0,oo). 

c) // {a,d) ^ {a^ — ad — 1 = 0}, then Spec{A4D) is infinite but it 
only intersects the real line in a finite number of points. 

Proof . By virtue of v), A4 is diagonalizable. We assume a + d > 0, 
so that 

{6±} C {Re{z) > 0}. 

Let 



y := a/4 — (a — dy and ■j± = a — d ± iy. 

Then 

^-(^■^)(^ ")(4 f 

Let := arg7+ and a := 1/ so that a = 1/ ^Jbl. Then 

EV{x) 27+ . , , ,_ , . Ill 1-^ s 

^ ^ — = [1 — cos(aa;j cos(a;a;)J — — ^ \ \ T~ I sm(axj sm(aa;j 

27+ /_ \n / 7^ a «\ . / N . /_ \ 

= 1 — cos(ax) cos(ax)J — — ^ — h — sm^axj sm(axj 

7+ V7+ a ay 

= 2e*^''[l — cos(ax) cos(ax)] — ( q^^'^— + ^ ) sin(ax) sin(ax) 

\ a a J 

= 2e*^''[l — cos(«x) cos(ax) — csin(ax) sin(ax)] 
= 2e'2^F(x), 

(21) 
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where F{x) and 

_ (ae'^^/a) + {ae-'^^ /a) 
2 

gi(2i9-2arg(a)) _^ g-i(2i?-2 arg (a)) 

~ 2 

= cos(2'(9 — 2arg(a)) = cos(2'(9 + arg6+) 

are as in lemma 

Let us show a). The hypothesis a — d < ensures — 1 < c < 1. 
Furthermore c = 1, if and only if 

Im(7^) _ Im(6+) 

Re(7i) ^"Re(6+)- 

The latter occurs, if and only if 

y{a -d) y 



(a — — 2 a + d 

By simplifying this identity, we gather that c = 1 for — ad — 1 = 
which is precisely our assumption. Then, lemma ^3|-b) and (pl|) 
complete the proof of a). 

For b), notice that since a — d > 0, the constant c is now such that 
— 1 < c < 1 and c = —1, if and only if 

Im(7^) Im(6+) 

Re(7i) ^"Re(6+)- 

Therefore a similar argument as for a) and lemma ^Sj-a) show this case. 
In order to prove c) use the fact that — 1 < c < 1 in 

i?5 \ {a^ - ac/ - 1 = 0}, 
lemma ^3|-c) and ■ 

Theorem 25. In the regions and Ri^, Spec(A4D) is infinite, and 
Spec {A^D) C {(r + iyof : r G M} + [0, oo) (a, d) G R^ 

Svec{AiD) C {-(r + i?/o)^ : r G M} + (-00,0] {a,d) e R^, 

where in both cases the constant yo > only depends upon {a,d). 

Proof . We show the result only for R^. According to iii), in this case 
< 6_ < 6+ and is diagonalizable. By expressing the trigonometric 
functions in exponential form, 

EV{x) = ki — ki cos(ax) cos(/?x) — k2 sin(Q;x) sin(/5x) 
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where ki,k2 G M and < f3 < a are constants we do not need to 
specify here. A similar argument involving Hadamard's theorem as in 



the proof of lemma ^ shows that Spec{A4D) is infinite. 

By putting x = r + iy where r, y G M, 7 := a + /5 > and 6 := 
a- /3 > 0, 

EV{r + iy) = ki + ^Q-^y^hr ^ ^^y^-^r^ ^ 

_ ^2 + ^1 Te'^^e*'^'' + e'^^e"*'^''! 
4 L J • 

Since 7 > 5 > 0, if we chose y » 0, the term e^^ dominates the 
expression and so \EV{r + iy)\ > c > for a suitable c independent of 
r. If we chose y << 0, the term e""^^ is the one that dominates and 
again \EV{r + iy)\ is large. This shows that all the zeros of EV{x) 
must be contained in a band {—yo < Im(a;) < yo}. ■ 

The above theorem does not rule out the possibility of negative eigen- 
values when ad < 0. We will see in the numerical examples, evidence 
of points in this region such that has indeed negative spectrum. 

With regard to finding the minimal yo. We will see in section |^ an 
argument involving Chebyshev polynomial that allows us to compute 
in closed form SpeclA^D) for a certain dense subset of R^. We will 
also illustrate this technique in various examples where the parabolic 
region is found explicitly. 

If {a,d) G Ri, the matrix A4 is not diagonalizable and so EV{x) 
is given by theorem |16| instead of theorem Nevertheless, similar 



techniques to the ones we have seen so far apply to this case. 
Lemma 26. Let 7^ c G M and let 

F{x) = + c[sin(x)]^ x G C. 

Then F{x) has an infinite number of zeros in the complex plane but 
only a finite number of them lie on M and on iM. 



Proof . See the proofs of lemma and theorem 25 



Theorem 27. Lei {a,d) G Ri. If {a,d) = (±1/2,^3/2), then 
Spec{A4^D) = {0}. Otherwise SpeclA^D) is infinite but it only in- 
tersects the real line in a finite number of points. 

Proof . lia-d = 2, 

+ 1 -1 ^ A A fb+ 0\ (I -1 
1 6+ - 1 i lo 1 A 1 h+ 1 1 



^4 

and ii a — d = —2 



1 

6+ - 1 -1 \ A -l\ fh+ 0\ /I 1 



1 h+ + l} VO 1 M 1 6+/ VO 1 
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Then 

o / , N 2 r / 

X 



X^ 



46i V 26 



sm 



a-d = ±2. 



The first statement follows from the fact that if (a, d) = (±1/2, =1=3/2), 
then 6+ = =Fl/2 and so the trigonometric term disappear. The second 



follows from lemma 26 



Notice that the curve — ad — \ = meets the region Ri at 
(±1/2, =p3/2). These are the only points where SpeclA^D) is finite. 
Since all self-adjoint operators with compact resolvent must have an 
infinite number of eigenvalues, is not similar to self-adjoint. All 
this suggests that for (a, d) in a small neighbourhood of these points, 
Spec{A4^D) must be highly unstable. In the next section we explore 
closely this idea. 

7. Some numerical results 

As mentioned previously, this section is devoted to investigating some 
aspects of the global spectral evolution of AD when we move the entries 
of the matrix A. To be more precise, we consider A = A4 (see section 
^) and compute Speci^A^D) as {a,d) moves along various lines inside 
i?i U i?3 U -R5 C M?. We also introduce a technique that allows us to 
find explicitly Spec {A4D) when (a, d) are in a certain dense subset of 
i?3 by computing the roots of certain polynomial G{w). 

Our first task is to decompose R3 into a disjoint union of curves in 
order to find the dense subset. For a > 1, let 

, , , a(a^ ± 1) ± v/(«^ - l)2a2 ± 4a2(a2 ± 1)2 
d^a) := 

and let 

A±(a) := {{a,d±{a)) : a > ±1}. 

Then 

i?3 n {a - < 0} = y A+(a) and R3 n {a - d > 0} = (J A-{a) 

The motivation for this decomposition is found by observing that for 

-1 
1 dJa) 



A, 



\/b+/b_ = a, where < 6_ < 6+ are the eigenvalues of A^. That is, 
A± are level curves of y^b^Jbl in the (a, (i)-plane. Notice that 

i?3 = IJ A+(a)uA_(a). 

The key idea behind finding G{w) is that for {a,d) G A±{a) where 
1 < a G Q, the zeros of the transcendental function are periodic in the 
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horizontal direction. We show how to construct this polynomiaL The 
transcendental function for A4^D is 



EV{x) = ki 



X 



cos — = COS 



+ 



— ^2 sin — = 1 sin 



X 



where ki and /c2 are two real constants depending upon a and d which 
we do not need to specify here. Since 

VK/b- = a=p/q, p,gGZ+, 

^/b± are rationally related and so the zeros of EV{x) appear periodi- 
cally in lines parallel to the real axis. By putting z = x/{q^/b^), 

EV{x) = ki[l — cos(pz) cos(g2;)] — k2 sin(pz) sin(g2;) 

= ki + ^L—h. cos[(p + q)z\ - cos[(p - q)z], 

where p — g<p + gGZ+. Standard computations show that, 

cos(mz) = Tm{cos{z)) m = 1, 2, . . . 

where a polynomial of degree m (the m*'^ Chebyshev polynomial of 
first order). Then by letting 

G{w) := fci + ^ ^^ T^p+q){w) - ^^^^T^p_q){w), 

EV{x) = 0, if and only if G{cos{z)) = 0. Hence all the zeros of EV{x) 
are of the form 



(±arccos(wo) + 2mT)qybj^ G C, n G Z 

where wq is a root of G{w). In this manner, Spec(yl4D) is generated 
by translations of the roots of G{w). 

Although the above method computes Spec(y44Z)) explicitly for 
(a, c?) G A-|-(a),l < a G Q, its numerical implementation for large 
p + g (> 20 in a PC) is highly unstable due to the well known instabil- 
ity of the roots of polynomials of high degree. Nevertheless, no other 
procedure tried so far, has proven to be more efficient for estimating 
large eigenvalues in i?3. Figures ||, ^ and H below were produced via 
this approach. 

7.1. Spectral behaviour of for (a, rf) close to ( — 1/2,3/2). 
By virtue of theorem 13, Spec(v44D) = {0} for (a, d) = (-1/2,3/2). In 
any small neighbourhood of this point, the spectrum of A^^D is infinite 
so high instability is to be expected. Since A^^D is holomorphic in a 
and d, every non-zero eigenvalue of A^^D either concentrates at zero or 
diverges to oo for (a, rf) (—1/2,3/2). We explore this phenomenon 
in some detail. 
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Figure 2. Evolution of the first 16 eigenvalues of A^D 
for a = —1/2 fixed and different values of d close to 3/2. 



According to theorem ^4|-a), if {a,d) G -R5 satisfy a — ad — 1 = 
and —2 < a — d < 0, 



Spec {A4D) = ^^-k\y[lm{b; 



1/2ni2 



where b+ as in section |6.3| . By taking a — — 1/2 and d 3/2, 

a + d+ J(a - d)2 - 4 
b+ = ^ ^ > 1/2 G M, 

— 1/2 

SO that Im(6_,_ ) ^ 0. Hence, all non-zero eigenvalues of A4D remain 
negative and escape to —00 as (a, d) G R5 approach the critical point 
on the curve — ad — 1 = 0. 

In general, not every eigenvalue of A4D need to be in the left hand 
plane when {a,d) is close to (—1/2,3/2). In figure |^ we consider the 
evolution of the first 16 eigenvalues of A^D for a = —1/2 fixed and 6 
different values of d from d = 1.6075 to 1.5035. The awkward choice 
of d correspond to the sensible values of a G Q; each pair (— 1/2, G 
A+(a) for a = 2,8/5,4/3,5/4,7/6,9/8. Notice that for large p,q the 
polynomial G{w) has p + q roots and nonetheless all these roots but 
lie on the same curve. This curve moves away from the origin and there 
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Figure 3. Evolution of the negative eigenvalue for a = 
— 1/2 and 100 different values of d linearly distributed on 
the segment [1.5012,1.6200]. 

is always a negative eigenvalue. The positive eigenvalues also escape 
rapidly to +oo and there are infinitely many of them. 

In figure |^ we isolate the negative eigenvalue for a = —1/2 against 
100 different values of d close to c? = 3/2. This provides indication of 
how rapidly it escapes to — oo. In order to produce this picture, we 
made use of the algorithm that Matlab provides to find the zero of 
EV{x) for X on the imaginary axis. Comparing with the comment we 
made earlier in section |673| , this provides points in such that A^D 
has a negative eigenvalue of arbitrarily large modulus. 

7.2. Non-real eigenvalues in Ri. 

We now explore the transition from real to non-real spectrum by con- 
sidering the spectral evolution of on the line 

{(0,rf) eRs : d>2} 

close to (0,2) E Ri. In figure || we show the first 23 eigenvalues of 
for a = and d = 2. We produced this graphic by reducing the 
equation EV{x) = to a single real variable and then making use of the 
algorithm that Maple provides to find zeros of real functions. According 
to theorem |2^, we know that Spec {A4D) is infinite but there is only 
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Figure 4. First 23 eigenvalues of A^D for a = and d = 2. 

finite intersection with the real line. As the picture suggests, in this 
case the origin seems to be the only real eigenvalue. 

Figure |^ shows the evolution of the first 17 eigenvalues (counting 
multiplicity) of A4 when a = for various different values of d from 
d = 3.3333 to 2.0139. Each pair {0,d) G A+(a) respectively for a = 
3, 5/2, 9/4, 2, 9/5, 3/2, 5/4, 9/8. The numerical evidence suggests 
that for d = 3.3333 the spectrum is close to the real line and each 
eigenvalue is of multiplicity 2. Each of these operators has infinitely 
many real eigenvalues. Unfortunately the method we employed to find 
the roots of G{w), is unable to deal with a finer partition of the d- 
interval. Nonetheless, the global behaviour of the spectrum can be 
appreciated, as d approaches to 2, each real eigenvalue eventually splits 
into two conjugate non-real single eigenvalues stabilizing close to the 
region in figure ^ (see the step d = 2.0139). Notice that there is no 
spectrum in the left hand plane and compare with theorem E2. 



7.3. Spectral evolution close to Rq. 

Another type of pecuhar behaviour can be observed as (a, d) G R3 ap- 
proach the region Rq, where the matrix A4 is singular and Spec {A4D) = 
C. Here we concentrate on the point (—1, 1) G -Re- 
Figure ^ shows the evolution of the first 100 eigenvalues of A^D 
(represented by dots) as {a,d) G A_,.(2) approaches to (—1,1) G -Re- 
Alongside we also picture the remaining eigenvalues (represented by 
crosses) that lie on the box [0, 2000] x [—300, 300] . A very similar 
behaviour occurs for {a,d) G A-i-(a) as {a,d) — > (=1=1, ±1) G -Re for 
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other values of a e Q. It can not be appreciated in the graph pro- 
vided but there are two conjugate eigenvalues whose real part is nega- 
tive. These eigenvalues approach to the origin as {a,d) (—1, 1). All 
the remaining spectrum concentrates on the real line suggesting that 
Spec (744/)) [0,00) as (a, d) — > (—1,1) this is in contrast with the 
fact that Spec (AD) = C at (-1, 1). 

Here we have chosen p = 2 and q = 1- This means that G{w) is only 
of order 3 and so the spectrum is always generated by 3 points. It is 
not difficult to show analytically that all three roots converge to and 
then rigorously prove that Spec{A4D) [0, 00). 
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Figure 5. Evolution of the first 17 eigenvalues (count- 
ing multipficity) of A4D for a = and d > 2 close to 
d = 2. 
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Figure 6. Evolution of the first 100 eigenvalues of A4D 
as (a, d) — > (—1, 1) G Bs on A+(2). The dots are the first 
100 eigenvalues while the crosses the remaining spec- 
trum. 
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